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ABSTRACT 

The dynamical behavior of a cantilever beam subject to parametric excitation is investigated. The method of 
multiple scale perturbation technique is applied to obtain the solution up to the second order approximation. At first we 
studied the system behavior without control and the effect of parameters at this case. The active control of negative 
feedback with linear and cubic velocity is used to vibration control. It can be seen that the linear control is more effective 
than the cubic control. All resonance cases are studied numerically to obtain the worst case. Stability of the steady state 
solution for the selected resonance case is studied applying Rung-Kutta fourth method and frequency response equation 
via Matlab 7. 0 and Maple 1 6. 
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1. INTRODUCTION 

Over the last few years, numerous papers have been dedicated to the control of resonantly forced systems 
in various engineering fields. Comparing active control to passive control, different studies prove that the later is 
more sophisticated than the former in various points. As in passive vibration absorbers, a physical device is 
connected to the primary structure, on the contrary, in active absorbers, the device is replaced by a control system 
of sensors, actuators and filters. 

Oueni et al. [1] studied a non-linear active vibration absorber coupled with the plant through user-defined 
cubic nonlinearities. Tondl and Nbergo [2] studied a dynamic absorber, which can move in the transverse or 
longitudinal direction, which attached to an externally excited simple pendulum. Nayfeh et al. [3] discussed 
dynamics of machining using quadratic and cubic stiffness of machine tools, which accounts for the regenerative 
effects. 

Hamdi and Belhaq [4] investigated nontrivial solutions and bistability in the Duffing oscillator with a 
delayed displacement feedback via the perturbation method. Eissa and co workers [5-7], studied the vibration 
reduction of a ship pitch-roll motion described by non-linear spring pendulum system under multi external, 
parametric excitations or both of them can be reduced using a longitudinal absorber. Eissa et al. [8-9] studied the 
active control of an aircraft tail subjected to multiharmonic excitation or multi-parametric excitation forces. They 
added many active controllers for this system. The system has controlled by negative (linear, quadratic) velocity 
feedback. Best active control of the system has been achieved via negative velocity feedback. 
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Y. S. Hamed et al. [10-12] studied the USM model subject to external and tuned excitation forces. The model 
consists of multi-degree-of-freedom system consisting of the tool holder and three absorbers (tools) simulating ultrasonic 
machining process. The advantages of using multi-tools are to machine different materials and different shapes at the same 
time. This leads to time saving and higher machining efficiency. Warminski et al. [13] analyzed the influence of the 
controller parameters, and also the hard or soft geometrical plant nonlinearities, to the dynamical response of the system 
undergoing NSC. Moreover, the effectiveness of the selected control algorithms to vibration suppression of a strongly 
nonlinear beam structure was compared and discussed. 

Lei et al. [14] applied an active control technique to coordinate a kind of two parametrically excited chaotic 
system. Oueini and Nayfeh [15] modelled the dynamics of the first mode of a cantilever beam with a second order 
nonlinear ordinary differential equation subjected to a principal parametric excitation, and a control law based on cubic 
velocity feedback is introduced. Pai and Schulz [16] studied the control of the first mode vibration of a stainless steel beam 
through negative velocity feedback to the dynamic system. 

Sayed and Kamel [17-18] studied the effect of different controllers on the vibrating system and the saturation 
control of a linear absorber to reduce vibrations due to rotor blade flapping motion. The stability of the numerical solution 
is investigated using both phase plane methods and frequency response equations. J. Warminski [19-20] studied apart from 
the regular dynamics, interactions between self- parametric or external excitations may lead to chaotic or hyperchaotic 
dynamics. 

Saeed et al [21] studied nonlinear time delay saturation-based controller for suppression of nonlinear beam 
vibrations. Gao and Chen [22] studied Active vibration control for a Bilinear system with nonlinear velocity time-delayed 
feedback. Yingli et al [23] studied Dynamic effects of delayed feedback control on nonlinear vibration floating raft 
systems. Gao, and Chen [24] studied Nonlinear analysis, design and vibration isolation for a bilinear system with 
time-delayed cubic velocity feedback. Amer and Ahmed [25] studied Vibration control of a nonlinear dynamical system 
with time varying stiffness subjected to multi external forces 

The objective of this work is to study dynamical behavior of a cantilever beam subject to parametric excitation 
under state feedback active control. The method of multiple scale perturbation technique is applied to obtain the solution 
up to the second order approximation. 

The equation of the dynamical behavior of a cantilever beam subject to parametric excitation is given by: 

X + sjuX + X + saX 3 +S/5XX 2 +syXX 2 = T +sX f cos($t) (1) 

where a dot denotes differentiation with respect to time^the damping coefficient, a the curvature nonlinearity 
coefficients, P and y are the inertia nonlinearity coefficients, / is the amplitude of the parametric excitation, (j) is the 

corresponding frequency, T = —£vX r , V — 1, 3. is the control V is the feedback gain. 

2. MATHEMATICAL TECHNIQUES 

The multiple scale perturbation method is conducted to obtain the approximate solution of Eq. (1). Assuming the 
solution in the form: 
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X(t;s)—x 0 (T 0 ,T i ) + £ x l (T 0 ,T x ) + s 1 x 2 + 0(s 3 ) 


( 2 ) 


where T 0 —t and T n — s n t , (n=l,2) are the fast and slow time scales, respectively. In terms of T () , T x and T 2 , 
the times derivatives transform are: 


d 

dt 


d 2 


-D 0 + sD ] ,—j-D 0 +2sD 0 D x +s (D l +2D 0 D 2 ) 
dt 


(3) 


Where D =— — , (n=0,l,3). 

dT, 


Substituting Eqs. (2) and (3) into Eq. (1) and equating the coefficients of the same power of e in both sides, we 


obtain: 


<9(£°):(Z) 0 2 +1)x 0 =0 (4) 

0(s) : ( D l + l)xj = -2D 0 D l x 0 - juD 0 x 0 - ax - /3 xqD 0 x 0 -yxlD 2 x 0 -vD 0 x 0 + f x 0 cos (</>t) (5) 

0 (^ 2 ) : 

2 2 2 2 2 

( D 0 + l)x 2 = -2 D 0 D x x x -D x x q - 2 D 0 D 2 x 0 - juD {) x x - /liD x Xq - 3ax 0 x x - /?x 0 £) 0 x 1 - Px q D x Xq 

- 2 y 0x o x 1J D o x o - yx^D\x x - 2 yxlD 0 D x x 0 - 2 /x 0 x x DqX 0 - 2 vD x x 0 + / x, cos(^) ® 

The general solution of Eq. (4) is 

x 0 (T 0 J x ) = A(T x )e iT ^A(T x )e~^ (7) 


where A is a complex function in T x , which can be determined from eliminating the secular terms at the next 
approximation. Substituting Eq. (7) into Eq. (5), then we have that: 


(D 2 + l)x, = (-2iD l A-ijuA-3aA 2 A-ij3A 2 A + 3yA 2 A-ivA)e iT ° + 

(-a A 3 -ip A 2 +yA 3 )e 3iT ° +^-(Ae i(1+m + Ae i( - >+m ) + cc 


( 8 ) 


where cc denotes the complex conjugate, after elimination the secular terms then x x can be written as: 
x I (7’ 0 ,7’ 1 ) = E x e 3iT ° + E 2 e i(X+ * )T ° + E 2 e i(A+ * )T ° + cc 
Substituting Eq. (7) and (9) into Eq. (6), then we have that: 


(9) 
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(Z) 0 2 +l)x 2 = [-D 2 A- 2/D, A-/uD x A-3aE x A 2 A l -ipE ] A 2 A' -p(D t A)A 2 -2[](D ] A)AA+\ lyE t A 3 A 
-4iy(D ] A)AA + 2iy(D ] A)A 2 -v(D l A)~E 2 fA]e iT(> +[-6iE l A 3 -3iE x juA 3 -6E t aA 4 A 

- 6iE, A 4 A+ 22 EjA 4 A - 2iy(D t A) A 2 - p{D,A)A 2 -3iE\vA 3 ]^ [-3EfiiA 5 -5iE { pA 5 
+1 \EjA 5 }e 5,T ‘ } +[-2iE 2 (\+</>)A-iE 2 ju(\+</>)A+6E 2 aA 2 A+4iE l p(\+</>)A 2 A 

- 2E 2 y(\ +</>f A 2 A - 2E 2 yA 2 A +iE 2 v(l+</>)]e i(1+m +[3E 2 aA 2 A-2iE 2 p(l+<fr)A 2 A 

- 2iE 2 pA 2 A- E 2 y(\ +<f>f A 2 A- 2E 2 yA 2 A]e i ~' +m EJA 3 e UH)T,i +[3 E 2 aA 3 (10) 

+iE 2 p(] +<f>)A 3 +2 iE 2 A 3 - E 2 y(l+<f>) 2 A 3 -2 E 2 yA 3 +^E l fA 3 ]e {3+m ‘ 

T E 2 /Ae i(1+2m 


Elimination the secular terms of Eq (10), then the solution of X 2 can be written as: 

x 2 (T 0 ,T l ) = E 4 e 3iT ° + E 5 e 5iT ° + E 6 e Kl+m + E 7 e i(A+m + E s e i(3 ~* )T ° + E 9 e i(3+m 


+ E l0 e i(1+2m +cc 


( 11 ) 


where Ej , (j = 1 , 1 0) are complex functions of T { . 

3. STABILITY ANALYASIS 

We study the different resonance numerically to see the worst resonance, one of the worst cases has been chosen 
to study the system stability. The selected resonance case (j) — 2 . In this case we introduce the detuning parameter a 
according to 

d) - 2 + so 


( 12 ) 

Substituting Eq. (12) into Eq. (8), then eliminating the secular and small divisor terms from x ± , we get the 


following 


-2iD l A-ijuA-3aA 2 A-ipA 2 A + 3yA 2 A-ivA + ^-Ae iaT ' = 0 

To analyze the solution of equation (13) it is convenient to express A in the polar form: 

A(T l ) = ffl(r,)e W(ri) 


(13) 


(14) 


where a and 6 and are unknown real-valued functions. Inserting Eq. (14) into Eq. (13) and separating real and 
imaginary parts, we have 


3 f 

ad’ = -(a-y)a 3 -^j-acos<p 


(15) 
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(16) 


where (p = oT x —76 


For steady state solutions a' =<p' = 0 ,6' = — a obtained 


/ cos (p = — (a -y)a 2 -la 


(17) 


/ sin<£> = 2ju + — f5 + 2v 


(18) 


From equations (17-1 8) the frequency response equation is given by: 


(| (a - y)a 2 - 2cr) 2 + (2 ju + ^a 2 + 2v) 2 - f 2 = 0 

3.1 Linear Solutions 

Now, to study the stability of the linear solution of the obtained fixed let us consider A in the form: 


(19) 


A(T l ) = ^(p -iqY 


ST , 


( 20 ) 


where p and q are real values functions in T x and consider 28 = <J . 

Substituting from Eq. (20) into the linear parts of Eq. (13) and separating real and imaginary parts, the following 
system of equations are obtained: 


P' + + v)p + (S ~\f)q = 0 

2 4 


( 21 ) 


q'+i-S -^)P + yO + v)q = 0 


( 22 ) 


-^C“ + v)-A -(£“/) 
(<j+2/) -2(^+v)-^ 


= 0 


(23) 


where: 


^ 2 =-(^ + v )± 2 V 64 ( 5 2 -/ 2 


The linear solution is stable in this case if and only if the real parts of V > 0, p > 0 , and otherwise it is unstable. 
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3.2 Non-linear Solutions 

To determine the stability of the fixed points, one lets 

a = a l0 +a n and (p = (p m +(p n (24) 


where a l0 , (p l0 are solutions of Eqs. (17), (18) anda n , (p n are perturbations which are assumed to be small 
comparing toa 10 , $? 10 . substituting Eq.(24) into Eqs. (15) and (16) using Eqs. (17) and (18) and keeping only the linear 
terms in a n , (p n we obtain: 


, ,-,U 3/3 2 V f 

«n = (— -«io--Mi+(-«io c° s Pio)<Pii 


(25) 


9 f f 

(pu = (--(« -r)a 10 +- — cos^ 10 )a u +(- — sin <p l0 ) <p n 
8 4 a w 4 a w 


(26) 


_ v 

2 8 10 2 


/ 


*10 COS <?10 


9 f f 

--(a-y)a w + - — cos^ 10 - — sin^ 10 

8 4rz in 4 a in 


= 0 


(27) 


Consequently, a non-linear solution is stable if and only if the real parts of both Eigen values of the coefficient 
matrix (27) are less than zero. 

4. NUMERICAL RESULTS 

The dynamical system (1) is solved numerically by using fourth order Runge Kutta method at non resonance case 
without control as shown in figure 1 from this figure the system is stable with amplitude about 0.15 which considered as 
basic case. From figure 2 we can see that the system at resonance case (f) — 2 without control the amplitude is increased to 
about 0.85 (six times) compared with the basic case in figure 1. The effect of negative linear feedback velocity at 
resonance (/) — 2 is shown in figure 3 the amplitude is reduced to minimum value. In figure 4 the effect of negative cubic 
feedback velocity at resonance (f) — 2 make the system is stable with amplitude about 0.15, which mean that the negative 
linear feedback control is more effective than the negative cubic feedback control. Figure 5, shows the effects of different 
parameters, we can see that the amplitude is monotonic decreasing function of linear damping /i , the damping coefficient, 
the curvature nonlinearity coefficients a, the inertia nonlinearity coefficients P and y as shown in Figure 5a-5d respectively. 
But the amplitude of vibration is monotonic increasing function of the excitation force amplitude / as shown in Figure 5e. 
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Figure 1: System Behavior without Control 
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Figure 2: System Behavior without Control at Resonance Case (f> = 2 
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Figure 3: System Behavior with Negative Linear Feedback Velocity at Resonance (j) — 2 
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Figure 4: System Behavior with Negative Cubic Feedback Velocity at Resonance (j) — 2 



Figure 5: Effect of Parameters without Control 


4.1 Frequency Response Curves 

The frequency response equation (19) is a non-linear algebraic equation in the amplitude a. This equation is 
solved numerically and the results are shown in the Figure 6, which representing the vibration of the amplitude a against 
the detuning parameter o for the given values of the other parameters. The steady state amplitude is monotonic deceasing 
function of the gain V , the linear damping \i } the curvature nonlinearity coefficients a, the inertia nonlinearity coefficients 
(3 and y as shown in Figure 6a, 6d, 6eand 6f respectively. But the state amplitude is monotonic increasing function of the 
excitation force amplitude / as shown in Figure 6c, which are a good agreement with the numerical solution in Figure 5. 
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Figure 6: Frequency Response Curve for f=4, a=0.4, p=l, y=6, p=0.5, v=0.5, t=1 
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5. COMPARISION BETWEEN NUMERICAL SOLUTION AND APPROXIMATE SOLUTION 

Now, we get a good agreement of the approximate solution obtained from frequency response equation and the 
numerical solution obtained by Runge-Kutta-fourth-order method as shown in Figure 7. 



Figure 7: Comparison between the Numerical Solution ( ) and the Approximate Solution ( ) 

6. CONCLUSIONS 


The dynamical behavior of a cantilever beam subject to parametric excitation is investigated. The method of 
multiple scale perturbation technique is applied to obtain the solution up to the second order approximation. At first we 
studied the system behavior without control and the effect of parameters at this case. The active control of negative 
feedback with linear and cubic velocity is used to vibration, control. All resonance cases are studied numerically to obtain 
the worst case. Stability of the steady state solution for the selected resonance case is studied applying Rung-Kutta fourth 
method and frequency response equation. From the above study the following may be concluded: 

• The worst resonance case is the sub harmonic resonance <f> = 2 , which leads to the amplitudes of the main system 
are increased to about 0.85 (6 times) of the basic case. 

• The linear feedback velocity is better than the cubic feedback velocity to reduce the amplitude of vibration. 

• The steady state amplitude is monotonic deceasing function of the gain The steady state amplitude is monotonic 
deceasing function of the gain V , the linear damping ji, the curvature nonlinearity coefficients a, the inertia 
nonlinearity coefficients P and y. But the state amplitude is monotonic increasing function of the excitation force 
amplitude /. 

• The is a good agreement between the numerical solution and the approximate solution. 
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